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ABSTRACT 



The aim of the present study is to test wh ether the cold accretion of gas 
through a "cosmic filament" (IMaccio et al.ll2006l ) is a possible formation scenario 
for the polar disk galaxy NGC 4650A. If polar disks form from cold accretion of 
gas, the abundances of the HII regions may be similar to those of very late-type 
spiral galaxies, regardless of the presence of a bright central stellar spheroid, with 
total luminosity of few lO^L©. 

We use deep long slit spectra obtained with the F0RS2 spectrograph at the 
VLT in the optical and near-infrared wavelength ranges for the brightest Hll 
regions in the disk polar disk of NGC 4650A. The strongest emission lines ([Oil] 
Hp, [OIII], Ha) were used to derived oxygen abundances, metallicities and the 
global star formation rates for the disk. The deep spectra available allowed us 
to measure the Oxygen abundances {12 + log{0 / H)) using the Empirical method 
based on intensities of the strongest emission lines, and the Direct method, based 
on the determination of the electron temperature from the detection of weak 
auroral lines, as the [OIII] at 4363 A. 

The Oxygen abundance measured for the polar disk is then compared with 
those measured for different galaxy types of similar total luminosities, and then 
compared against the predictions of different polar ring formation scenarios. 
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The average metallicity values for the polar disk in NGC 4650A is Z = 0.2Zq, 
and it is lower that the values measured for ordinary spirals of similar luminosity. 
Moreover the gradient of the metallicity is flat along the polar disk major axis, 
which implies none or negligible metal enrichment from the stars in the older 
central spheroid. 

The low metallicity value in the polar disk NGC 4650A and the fiat metallicity 
gradient are both consistent with a later infall of metal-poor gas, as expected in 
the cold accretion processes. 

Subject headings: galaxies: individual (NGC4650A) - galaxies: formation - 
galaxies: evolution - galaxies: interactions - galaxies: peculiar - galaxies: abun- 
dances. 



1. Introduction 



The hierarchical, merger -dominated pictu re of galaxy formation is based on the Cold 
Dark Matter (CDM) model (jCole et al.l 120001 ). which predicts that the observed galaxies 
and their dark halo (DH) were forrned th rough a repeated merging process of small systems 
(jPe Lucia et al.l l2006t iGenel et al.l l2008l ). In this framework, major and minor mergers 
of disky systems do play a major role i n the formation of spheroid and elliptical galaxies 
(iNaab et al.ll2007l : iBournaud et al.l 120071). in all environments and from the Local Group to 
high-redshift universe (IConselice et al.ll2003l ). The gas fraction is a key parameter in the 
physics of such gravitat ional interactions: if it is high enough, an extended and m assive disk 
structure can survive (jSpringel fc HernquistI l2005l : iRobertson fc BuUockl |2008| ) . Galaxies 
can get their gas through several interacting processes, such as smooth accretion, stripping 
and accretion of primordial gas, which are equally important in the growth of galaxies. 
Recent theoretical works have argued that the accretion of external gas from the cosmic we b 
filaments, with inclined angular momentum (jPave et al.l 1200 ll . ISemelin fc Combed l2005l ) . 
might be the most realistic way by which gala xies get their gas. This process may also 



explain the build-up of high redshift disk galaxies (iKeres et al.ll2005l lKeresll2008l . iBrook et al. 



20081 . iDekel et al.ll2009l iBournaud Sz Elmegreen l2009h . The relative share of all gravitational 
interactions depends on the environments and it drives many morphological features observed 
in galaxies, such as bars and polar rings. 



"'^Based on the data acquired at the VLT with F0RS2, during the observing runs < 078. i? — 0580(A) > 
and < 079.B0177(A) >. 
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Galaxies with polar rings (PRGs) generally contain a central featureless stellar spheroid 
and an elongated structure, the "polar ring" (made up by gas, stars and du st), which orbits 



in a n early perpendicular plane to the equatorial one of the central galaxy ( iWhitmore et al. 
1990l ). The decoupling of the angular momentum of the polar structure and the central 



spheroid cannot be explained by the collapse of a single protogalactic cloud: thus a "sec- 
ond event" must have happened in the formation history of these systems. This is the 
reason why studying PRGs promise to yield detailed information about many of the pro 



cesses at work during galaxy i nteractions and merging (llodice et al. 



2002 



lodice et al. 2002 



Reshetnikov fc Combeslll994l . iReshetnikov et al.lbood ISournaud fc Combesll2003h 



The debate on the origin of PRGs is still open and two main processes have been 
proposed i) a. major dissipative merger or ii) gas accretion. In the merging scenario, the 
PRG result s from a "polar" merge r of two disk galaxies with unequal mass, ( iBekkil 119971 : 
Bekkil Il998t iBournaud et al.l 120051 ): the morphology of the merger remnants depends on 
the merging initial orbital parameters and the initial mass ratio of the two galaxies. In 
the accretion scenario, the polar ring may form by a) the disruption of a dwarf companion 
galaxy orbitating around an early-type system, or by b) the tidal accretion of gas stripping 
from a disk galaxy outskirts, c a ptured by an early-type ga l axy on a parabolic e ncounter 



(IReshetnikov fc Sotnikoval Il997 



Bournaud &: Combes 



20031 : iHancock et al.ll2009[ ). In the 



latter case, the total amount of accreted gas by the early-type object is about 10% of the gas 
in the disk donor galaxy, i.e. up to 10^ Mq. Both major merger and accretion scenarios are 
able to account for many observed PRGs morphologies and kin ematics, such as the ex istence 
of both wide and narrow rings, helical rings and double rings (IWhitmore et al.lll990l ). 



Very recently, a new mechanism has been proposed for the formation of wide disk- 
like polar rings: a long-lived polar structure may form through col d gas accretion alo ng a 
filament, extended for ~ 1 Mpc, into the virialized dark matter halo (IMaccio et al.ll2006l ). In 
this formation scenario, there is no limits to the mass of the accreted r naterial, thus a ver y 
massive polar disk may develop either around a stellar disk or a spheroid. iBrook et al.l (120081 ) . 
by using high-resolution cosmological simulations of galaxy formation, have confirmed and 
strengthened the formation scenario proposed by iMaccio et al.l (120061 ). In this case polar 
disk galaxies can be considered as extreme examples of angular momentum misalignment 
that occurs during the hierarchical structure formation. In the merging history, an inner disk 
formed first after the last major merger of two galaxies with a 1:1 mass ratio, then, due to 
the high gas fraction, the galaxy rapidly forms a new disk whose angular momentum depends 
on the merger orbital parameters. At later times, gas continues to be accreted along the disk 
which may be also in a plane perpendicular to the inner disk. The morphology and kinematics 
of one simula ted object, in both simulations, are similar to those observed for NGC4650A: 
in particular, iBrook et al.l ( 120081 ) found that such formation mechanism can self-consistently 
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explain both morphology and kinematics of central spheroid and polar structure, and all 
the observed features (like colors and colors gradient, longevity, spiral arms, HI content and 
distribution). 

NGC 4650A is the prototype for PRGs (see Fig. [1]). Its luminous components, inner 
spheroid and polar structure were studied in optical and near-infrar ed (NIR) photometry, 
spectroscopy, and in the radio emission, HI 21 cm line and continuum (lArnaboldi et al.lll997l . 
Gallagher et aPbood . Ilodice et al. l booi ISwaters fc Rubinlboosl . Ilodice et allbood ). 



The polar structure in NGC 4650A is a disk, very similar to that of late-type spirals or 
LSB galaxies, rather than a ring. The polar disk stars and dust can be reliably traced to 
~ 1.2 kpc radius from the galaxy nucleus, and the surfac e brightness profiles have an expo- 
nential decrease (Ilodice et al.ll2002UGallagher et al.ll2002l ). Furthermore, the rotation curves 
measured from the emission and absorption line optical spectr a are consistent with tho se of a 
disk in a differential rotation rather than that of a narrow ring (ISwaters fc Rubirul2003l ). This 



is also confirmed by the HI 21 cm observations (lArnaboldi et al.lll997l ) which show that the 



gas is five times more extended than the luminous polar structure, with a position-velocity 
diagram very similar to those observed for edge-on disks. The polar disk is very massive, 
since the total HI mass is about lO^^M©, which added to the mass of st ars, makes the mas s 
in the polar disk comparable with the total mass in the central spheroid (Ilodice et al.ll2002l ). 
The morphology of the central spheroid resembles that of a low-luminosity early-type galaxy: 
the surface brightness profile is described by an exponential law, with a small exponential 
nucleus; its inte grated optical vs. NIR colors are similar to those of an intermediate age stel- 
lar population ( Ilodice et al.ll2002l . iGallagher et al.ll2002l ). New high resolution spectroscopy 
in NIR on its photometric axes suggests that this component is a nearly-exponential oblate 
spheroid supported by rotation (Ilodice et al.ll2006l ). 



These new kinematic data, together with the previous studied, set important constraints 
on the possible formation mechanisms for NGC4650A. 

Because of the extended and massive polar disk, with strong Ha emissions, NGC 4650A 
is the ideal candidate to measure the chemical abundances, metallicity and star formation 
rates (SFR) via spectroscopic measurements of line emissions along the major axis of the 
polar disk. The goal is to compare the derived values for the metallicity and SFR with those 
predicted by the different formation scenarios. As we shall detail in the following sections, if 
the polar structure forms by accretion of primordial cold gas from cosm ic web filament, we 
expect the disk to have lower metallicities of the order of Z ~ I/IOZq (lAgertz et al.ll2009l ) 
with respect to those of same- luminosity spiral disks. 



We shall adopt for NGC4650A a distance of about 38 Mpc based on Hq = 75 km s ^ Mpc 
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and an heliocentric radial velocity V = 2880 km s ^, which implies that 1 arcsec = 0.18 kpc. 

2. Observations and data reduction 

Spectra were obtained with F0RS2@UT1, on the ESO VLT, in service mode, during 
two observing runs: 078.B-0580 (on January 2007) and 079.B-0177 (on April 2007). F0RS2 
was equipped with the MIT CCD 910, with an angular resolution of 0.25" pixel"^. The 
adopted slit is 1.6" wide and 6.8' long. Spectra were acquired along the North and South 
side of the polar disk, at P. A. = 152° (see Fig. [T]), in order to include the most luminous 
HII regions in polar disk. The total integration time for each direction is 3 hours, during the 
078.B-0580 run and 2.27 hours during 079.B-0177 run, respectively, with an average seeing 
of 1.2". 

At the systemic velocity of NGC 4650A, to cover the red-shifted emission lines of the 
[0/J]A3727, [i/^]A4340, [OJJJ]A4363, [0///]AA4959, 5007, [Hp]Xmi, [NII]X5755 the grism 
GRIS-600B+22 was used in the 3300 — 6210 A wavelength range, with a dispersion of 
50 A/ mm (0.75 A/pix). In the near-infrared 5600 — llOOOA wavelength range, the grism 
GRIS-200I+28 was used, with a dispersion of 162 A/mm (2.43 A/pix) to detect the fainter 
[5//]AA6717, 6731 and [5///]AA9068, 9532 emission lines, with a S/N > 20 and the brighter 
[i^Q,] A6563 emission line, with a. S/N > 150. In this wavelength range, where the sky back- 
ground is much more variable and higher with respect to the optical domain, to avoid the 
saturation of the sky lines a larger number of scientific frames with short integration time 
(850 sec) were acquired. 

The data reduction was carried out using the CCDRED package in the IRAfI {Image 
Reduction and Analysis Facility) environment. The main strategy adopted for each data- 
set included dark subtractioiifl, fiat-fielding correction, sky subtraction and rejection of bad 
pixels. Wavelength calibration was achieved by means of comparison spectra of Ne-Ar lamps 
acquired for each observing night, using the IRAF TWODSPEC.LONGSLIT package. The 
sky spectrum was extracted at the outer edges of the slit, for r > 40 arcsec from the 
galaxy center, where the surface brightness is fainter than 24'mag/ arcsec'^, and subtracted 
off each row of the two dimensional spectra by using the IRAF task BACKGROUND in 
the TWODSPEC.LONGSLIT package. On average, a sky subtraction better than 1% was 



^IRAF is distributed by the National Optical Astronomy Observatories, which is operated by the Asso- 
ciated Universities for Research in Astronomy, Inc. under cooperative agreement with the National Science 
Foundation. 

■^Bias frame is included in the Dark frame. 
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achieved. The sky-subtracted frames, both for North and South part of the polar disk, were 
co-added to a final median averaged 2D spectrum. 

The final step of the data-processing is the flux calibration of each 2D spectra, by using 
observations of the standard star LTT4816 and the standard tasks in IRAF (STANDARD, 
SENSFUNC and CALIBRATE). The flux calibration is very important in the present study 
because we need to "align" two spectra, which cover different wavelength range and taken 
in different times. Thus, we checked and obtained that the calibrations for the spectra in 
both spectral range were consistent. Fig. |2] shows the ID flux calibrated spectra of the 
spectrophotometric standard star used to calibrate the spectra in the whole range 3300 — 
11000 A. To perform the flux calibration we extracted a 1-D spectrum of the standard star to 
find the calibration function; then we extracted a set of 1-D spectra of the galaxy summing 
up a number of lines corresponding to the slit width. Since the slit width was 1.3" and the 
scale of the instrument was 0.25"/pix, we collapsed seven lines to obtain each 1-D spectrum. 
Finally we applied the flux calibration to this collection of spectra. 

Furthermore, we compared our flux calibrated spectra with oth ers acquired at the Siding 



Spring Observatory with the Double Beam Spectrograph (DBS) (IButtiglione et al.ll2006l ). 
The DBS has a dichroic that slits the light in a red and a blue arm, therefore the flux 
calibration with standard stars can be done simultaneously for the red and the blue arms. 
We used these spectra to check for any difference in the flux calibrations, finding that our 
flux calibrated spectra, both of the template star and of the galaxy, turn out to be consistent 
with them. 

The wavelength and flux-calibrated spectra are shown in Fig. [3] and |H In the blue 
spectrum (top panel) are clearly visible a number of emission lines: Hp, H^, [0//]A3727 
and [O J/J] A A4959, 5007, while in the red one (bottom panel) we have JJq, (blended with the 
[A^JJ]A6583 line appearing in the red wing of JT^), [^JJ]AA6717, 6731 and JJ J] AA9069, 9532. 
From a two-Gaussian fit to the combined emission, we estimate the line ratio [A^JJ] A6583/( JJq, -|- 
[A^JJ]A6548 ~ 0.1. For this reason the JJ^ flux is that measured as the total flux in the line 
reduced by the contribution of [A^JJ]A6583. The observed emission lines and their uncor- 
rected and reddening-corrected fluxes relative to JJ^ are listed in Tab. [2] and Tab. [31 

Since ground-based near-infrared spectroscopy is affected by the strong and variable ab- 
sorption features due to the Earth's atmosphere, we accounted for this effect in our spectra, 
in which the telluric absorption bands are clearly visible around 7200, 8200 and 9300 A. In 
order to perform this correction we used the "telluric standard star" LTT4816, near both 
in time and air mass to the object, and the IRAF task TELLURIC; by fitting the contin- 
uum, telluric calibration spectra are shifted and scaled to best divide out telluric features 
from data spectra. The ratios of the uncorrected fluxes relative to those obtained applying 
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the average telluric correction are the following: [S I I]\Q7VJ uncorrected/ [S I corrected = 

[SII]X6731 uncorrected/ [5 J J] A6 73 1 

corrected 

1.7 and [Sin]X9069uncorrected/[Sni]X9069 

corrected 

[S 1 1 1]X9532uncorrected/ [S 1 1 1]X9532corrected = 0.6. 



2.1. Measurement of emission-lines fluxes 

The fluxes of the above mentioned emission lines were measured using the IRAF SPLOT 
routine, that provides an interactive facility to display and analyze spectra. The if/j is 
evaluated for r > 10 arcsec, where only the emission line is present; for lower distances, 
i.e. where stars relative to the spheroid also contributes to the spectra, the Hp is also in 
absorption. We evaluated flux and equivalent width by marking two continuum points around 
the line to be measured. The linear continuum is subtracted and the flux is determined by 
simply integrating the line intensity over t he local fitted continuum. Th e errors on these 



quantities have been calculated, following iPerez-Montero fc Diazl ( l2003l ). by the relation 



(Ti = (7cA^^/^[l + EW/ (NA)]^^"^ , were ai is the error in the line flux, cJc is the standard 
deviation in a box near the measured line and represents the error in the continuum definition, 
N is the number of pixels used to measure the flux, EW is the equivalent width of the line 
and A is the wavelength dispersion in A/pixel. The errors relative to each emission line 
fluxes are listed in Tab. [2] and Tab. [31 



2.2. Reddening correction 



Reduced and flux calibrated spectra and the measured emission line intensities were cor- 
rected for the reddening, which account both for that intrinsic to the source and to the Milky 
Way. By comparing the intrinsic Balmer decrements, Ha/Hj^ = 2.89 and H^/H/s = 0.468, 
predicted for large optical depth (case B) and a temperature of lO'' K, with the observed 
one, we derived the visual ext inction A(V) a nd th e color excess E{B — V), by adopting 
the mean extinction curve by ICardelli et al.l (119891 ) A{X)/A(y) = a{x) + b{x)Rv, where 
Rv[= A(y)/E{B — V)] = 3.1 and x = 1/X. In order to estimate reddening correction for all 
the obse rved emission l ines, w e used the complete extinction curve into three wavelengths 
regions (ICardelli et al.l (119891 )): infrared (A > 0.9/im), optical/NIR (0.9/im > A > O.Sfim) 
and ultraviolet (0.125/im > A > O.lOyum), which are characterized by different relations of 
a(x) and b(x). All the emission lines in our spectra are in the optical/NIR range, except 
for the [5'//J]A9532, that falls in the infrared range, so we used the average i?y-dependent 
extinction law derived for these intervals to perform the reddening correction. 
We measured the fluxes of if^, and lines at each distance from the galaxy center and 
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for each spectra (North sht and South sht), than we derived the average observed Balmer 
decrements, which are the following: 
Ho,/ Hp = 2.40 ± 0.01 
H^/H(s = 0.41 ± 0.01 

while the color excess obtained by using these observed decrements are: 
[E{B-V)]h^/h, = 0.20 ± 0.004 
[E{B - V)]h,/h, = 0.25 ± 0.012. 

Such values of E(B-V) are used to derive the extinction Ax, through the Cardelli's law; in 
particular, the [E{B — V)]h^/His and [E{B — V^)]hq/h^ are used respectively for the red- 
dening correction in optical and NIR wavelength range. The corrected fluxes are given by 
F^JF^^^ = F„4^/Fff 10°-^'^^"^^/3'. Observed and reddening-corrected emission line fluxes 
are reported in Tab. [2] and Tab. El 



3. Oxygen abundances determination 



The main aim of the present work is to derive the chemical abundances in the polar 
disk of NG C4650A: in what fo ll ow we evaluate the oxygen abu ndai ices following th e methods 



outlined in IPagel et al.l ( 119791 ). [Diaz fc Perez-Monterd ( 20001) and 



Osterbrock 



ilvugini (I2OOII). r e 



(jl989|) andlAUenl (119841 ) 



errin g 



to them as Empirical methods, and those introduced by 
or Direct methods. 

The large dataset available (presented in Sec. [2]) let us to investigate both the Empirical 
methods, based on the intensities of easily observable lines (Sec. 13. ip . and the Direct method, 
based on the determination of the electron temperature, by measuring the intensities of the 
weak auroral lines (Sec. 13. 2p . 

As described in details in the next sections, we have derived the oxygen abundance parameter 
12 -|- log{0 / H) along the polar disk, by using both the Empirical and Direct methods. 



3.1. Empirical oxygen and sulphur abundance measurements 



The Empirical methods are based on the cooling properties of ionized nebulae which 
translate into a relation between emission-line intensities and oxygen abundance. Several 
abundance calibrators have been proposed based 011 differe nt emission-line ratios: amon g the 
other, in this work we used R2?, (jPagel et al.l 119791 ). 5*23 (IDiaz fc Perez-Monterd I2OOOI ) and 
the P-method (|Pilvuginll20011). The advantages of different calibrators have been discussed by 
several authors (jPfaz fc Perez-Monterdl2000l . iKobulnicky fc Zaritskylll999l . iKobulnicky &: Kewley 
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2004 IPerez-Montero fc Dfajboosl . iKewlev fc Ellisonlboosl iHidako-Gamez fc Ramfrez-Fuentes 



20091 ) ■ 



The method proposed by iPagel et ah Jl979h is based on the variation of the strone; 
oxygen hnes with the oxygen abundance. iPagel et ahl (119791 ) defined the "oxygen abundance 
parameter" R23 = ([0//]A3727 + [0///]AA4959 + 5007) /H^, which increase with the oxy- 
gen abundance for values lower than 20% of the solar one and then reverses his behavior, 
decreasing with increasing abundance, due to the efficiency of oxygen as a cooling agent 
that lead to a decreasing in the strength of the oxygen emission lines at higher metallicities. 
Three different regions can b e identified in the trend of R23 with the oxygen abundance 
(jPerez-Montero fc Dfaal2005l ): a lower branch [12 + log{0 / H) < 8.1), in which i?23 increase 
with increasing abundance, an upper branch (12 + log{0/H) > 8.4), in which the trend 
is opposite, and a turnover region (8.1 < 12 + log{0 / H) < 8.4). While the upper and 
lower branches can be fitted by regression lines, the turnover region seems to be the most 
troublesome, in fact object with the same value of -R23 can have very different oxygen abun- 
dances. The R23 parameter is affected by two main problems: the two-valued nature of the 
calibration and the dependence on the degree of ionization of the nebula. 

To break the dege n eracy that affects the metallicity for values of 12 + log{0/H) > 8.0, 
Diaz fc Perez-Monterd (I2OOOI ) proposed an alternative calibrator, based on the intensity of 
sulphur lines: 523 = ([5//]AA6717 + 6731 + [5///]AA9069 + 9532)/iJ;3. The [SII] and [SIII] 
lines are analogous to the optical oxygen lines [Oil] and [OIII], but the relation between 
the " Sulphur abundai ice parameter 5*23" and the o xyge n abundance remain single va lued up 
to solar metallicities ( Diaz fc Perez-Montero 200ol and Perez-Montero fc Diaz 2005 ). More- 
over, because of their longer wavelength, sulphur lines are less sensitive to the effective 
temperature and ionization parameter of the nebula, and almost independent of redden- 
ing, since the [SII] and [SIII] lines can be measured nearby hydrogen recombination lines. 
Diaz fc Perez-Monterd (120001 ) have attempted a calibration of oxygen abundance through 
the sulphur abundance parameter, obtaining the following empirical relation: 



12 + log{0/H) = 1.53logS23 + 8.27 (1) 

For the polar disk of NGC4650A, by using the emission line fluxes given in Tab. [2] and 
Tab. [31 we derived both R23 and S'23 parameters, which are listed in Tab. HI Since there are 
few measurements available for the Sulphur lines, we have obtained an average value of S'23 for 
the polar disk, by using the following telluric corrected fluxes ratios: [S 1 1]\6717 / Hjs = 0.18, 
[SII]X6731/H^ = 0.2, [SIII]\9069/Hfi = 0.43 and [SIII]X9532/Hf} = 0.25, obtaining 
logS23 = 0.024 ± 0.030. The average oxygen abundance, derived by adopting the Eq. [H 
is the following: 12 + log{0/H)s = 8.3 ± 0.2. 
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Pilyugid ( 120011 ) realized that for fixed oxygen abundances the value of X23 = logR23 
varies with the excitation parameter P = R3/R23, where -R3 = 0///[4959 + 5007]/if^, and 
proposed that this latter parameter could be used in the oxygen abundance determination. 
This method, called " P-method" , propose to use a more general relation of the type 0/H = 
f{P, R23), compared with the relation 0/H = /(-R23) used in the R23 method. The equation 
related to this method is the following 

10 , i?23 + 54.2 + 59.45P + 7.31P2 

12 + log{0/H)p = eo7 + 6.71P + 0.371P2 + o.243i?23 

where P = R3/R23. It can be used for oxygen abundance determination in moderat ely high- 



metallicity HII regions with undetectable or weak temperature-sensitive line ratios ( iPilyugin 



2OOII) 



We used also the " P-method" to derive the oxygen abundance in the polar disk of NGC4650A: 
the values of 12 -|- log{0/H)p at each distance from the galaxy center are listed in Tab. HI 
and shown in Fig. [5] and Fig. |6l The average value of the oxygen abundance is 12 -|- 
log{0/H)p = 8.2 ± 0.2, which is consistent with the values derived by using the Sulphur 
abundance parameter (Eq. [1]). 

The metallicities corresponding to each value of oxygen abundances giv en before have 



been estimated. We adopted 12 + log{0 / H)q = 8.83 = Aq and Zq = 0.02 (lAsplund et al. 



20041). Given that Zngc4650A ~ KZq and K = loi^NGC4650A-AQ]^ obtain a metalhcity 
for the HII regions of the polar disk in NGC4650A Z 0.004 which correspond to 
Z ~ (0.2 ±0.002)^0. 



3.2. Direct oxygen abundance measurements 

The electron temperature is the fundamental parameter to directly derive the chemical 
abundances in the star-forming regions of galaxies. In a temperature-sensitive ion, with 
a well separated triplet of fine-structure terms, electron temperature and electron density 



Osterbrock 


1989: 


Allen 


1984) 



in detail below, the oxygen abundances are function of both the emission line fluxes and 
electron temperatures: this relation is the basis of the commonly known as direct method or 
Te method. 

Usually, it happens that not all these lines can be observed in the spectra, or they are affected 
by large errors. Thus, some assumption were proposed for the temperature structure through 
the star-forming region: for HII galaxies it is assumed a two-zone model with a low ionization 
zone, where the OH, Nil, Nell and SII lines are formed, and a high ionization zone where the 
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OIII, NIII, Nelll and SIII lines are formed (see, e.g. ICampbell et al.lll986l : lGarnettlll992h . 
Photoionization models are then used to relate the temperature of the low ionization zone ty 



to t.s, the temperature of the high ionizat i on zone (iPaee. 



20031 : IPilvurin et all 120061 : IPilvudnl 120071 : IPilvurin et al. 



et al JI1992I : iPerez-Montero fc Diaz 
2009h . 



For NGC4650A, we aim to derive the oxygen abundance of the polar disk directly by the 
estir nate of the 0^ ~ ^ and ions electron temperatures. According to llzotov et al.l (120051 ) 
and iPilyugin et al.l (120061 ) , we have adopted the following equations for the determination of 
the oxygen abundances, which are based on a five-level atom approximation: 



12 + log{OIII/H) = /o^([4959 + 5007]/Hp) + 6.2 + 1.251/^3 - O.SS/oc/tg - O.OMtg (3) 



12 + log{On/H) = log {[3727] /Hp) + 5.961 + 1.676/^2 - 0A0logt2 - 0.034t2 + log{l + 1.35x2) 

(4) 

where and t2 are the electron temperatures within the [OIII] and [Oil] zones respec- 
tively, in units of 10^ K; X2 = lO'^A^^gt^^'^, where Ne is the electron density. The total oxygen 
abundance is 0/H = OIII/H + OII/H. 

The electron temperatures and the electron density, needed to solve the above re- 
lations, have been derived using the task TENDEM of the STSDAS package in IRAF, 
which solves the equations of t he st a tistical equilibrium w ithin the five-level atom ap- 



proximation (jPe Robertis et al.l 119871 : IShaw fc Dufoun 119951 ). According to this model 



both ^3 and ^2 are function of the electron density and of the line ratios Rqiii = 
[OJ/J] A (4959 + 5007)/[OJJJ]A4363 and Ron = [OJJ]A3727/[OJJ]A7325, respectively. For 
NGC4650A, has been derived from the line ratio [5'//]A6717/[5J/]A6731, which has the 
average value of 0.83 ± 0.06, through the whole disk extension, for a default temperature of 
10^ K. This let to N^. ~ 1217 cm~^. Even though the auroral line [0///]A4363 is usually 
very weak, even in the metal-poor environments, the large collecting area of the 8m telescope 
and the high signal-to-noise spectra, let us to measure such emission line along the whole 
extension of the polar disk (see Tab. [2]), which let us to estimate Rom- 
The electron temperature = f {Rqiii, Ng) has been calculated by putting Ne and Rqiii as 
inputs in the TENDEM task: given the large spread of Rqih, we adopted the average value 
in the following three bins, < -R < 20 arcsec, 20 < i? < 40 arcsec, 40 < -R < 60 arcsec. The 
distribution of in each bin is shown in Fig. [71 the electron temperature is almost constant 
till about 40 arcsec and tends to increase at larger distances. In the same figure there are also 
plotted the values of ^2 (at the same distances) obtained by adopting the following empirical 
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t2 — ts relation, derived by iPilyugin et al.l (|2009[ ): t2 = 0.264 + 0.835^3. The empirical ^2 — is 
relation has been inve stigated and debated for a long time and several forms are given in 
the literature; that by iPilyugin et al.l (|2009[ ) is the most recent work on this subject, where 
they derived the t2 ~ relation in terms of the nebular = [O///] A(4959 + 5007) / Hj^ and 
i?2 = [O I I]X3727 / H p line fluxes. They found that such relation is valid for HII regions with 
a weak nebular R3 l ines (loqR^ > . 5) an d it turns to be consistent with those commonly 
used (see Fig. 14 in IPilyugin et al.l ( l2009l ) paper and reference therein). The polar disk of 
NGC4650A, since logR^ = 0.5 ± 0.2 (this is the average value on the whole disk extension), 
is located in the upper validity limit of the ^2 — ^3 relation. 

This is what is usually done since the auroral line [0//]A7325 is too weak to obtain a good 
enough estimate of the electron temperature ^2- For the polar disk of NGC4650A, the line 
fluxes of [O/J] A(7320 + 7330) are strong enough to be measured at about 30 arcsec and 
the average value of [O/J] A(7320 + 7330)/i7/3 = 0.04 ± 0.01. The average value of the ra- 
tio Ron = [0//]A3727/[0//]A7325 is given as input in the TENDEM task, and we obtain 
t2 = (1.092 ± 0.20) X 10^ K, which is added to Fig. [71 This value is lower than that de- 
rived by using photoionization models, even if the observed and theoretical values of ^2 are 
consistent within errors: such difference could be due both to the uncertainties in the line 
fluxes measurements, made on so few and weak lines, and to the validity limit of the ^2 — 
relation. 

To obtain the oxygen abundance for the polar disk, we decided to adopt both estimates for t2, 
by using the average value inside 40 arcsec for that derived by photoionization models, which 
is t^""^ = (1.7±0.2) X 10^ K, and, consistently, an average value for = (2.0±0.2) x 10^ K at 
the same distances. By adopting the Eq. [3] and Eq. Hj, the oxygen abundance 12 + log{0 / H)t 
for the polar disk is shown in Fig. [8] and listed in Tab. [51 By using the value of t2 de- 
rived by direct measurements of Oil line fluxes, the metallicity is higher with respect to 
that derived by using the value of ^2 from photoionization models: the average value is 
12 + log{0/H)T^^^ = 8.4 ± 0.1 and 12 + log{0 / H)t^^^ = 7.6 ± 0.5 for the two cases respec- 
tively. The large error, that are derived by propagating the emission line intensity errors 
listed in Tab. [2] and Tab. [3l which is about 0.5, let both estimates consistent. 
It is important to point out that the oxygen abundance derived by direct measurements of 
on line fluxes for t2 is much more similar to the value derived by empirical methods, which 
is 12 -|- log{0/H) = 8.3, than the value derived by using photoionization models for t2: this 
may suggest that NGC4650A could be out of the range of validity of these models. 
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4. Discussion: use of the chemical analysis to constrain the galaxy formation 

As discussed above, we derive the oxygen abundance, 12 + log{0 / H) to be 8.2 ± 0.1, 
by using both the empirical and direct methods. In the following sections we will discuss 
the results obtained by the present work, how they could reconcile with the predictions by 
theoretical models and, finally, we address the main conclusions of this study. 



4.1. Results 



Here we discuss how the average value of metallicity and its distribution along the polar 
disk compares with those typical for other late-type disk galaxies and PRGs. 

- Metallicity-Luminosity relation 
The mean values for the oxygen abundance along the polar disk, derived by the empirical 
(see Sec. 13. ip and direct methods fsee Sec. 13. 2p. are c ompared with those for a sample of 
late- type disk galaxies by iKobulnicky &: ZaritskyI (119991 ). as function of the total luminosity 
(see Fig. For all galaxies of the sample, and in particular for spiral galaxies, the total 
luminosity is that of the whole system, i.e. bulge plus disk, therefore also for NGC4650A 
we have accounted for the to tal luminosity of the ga l axy ( Mr = —19.3, evaluated by using 
the same value of Hq used by IKobulnicky fc ZaritskyI (119991 ) in order to compare NGC4650A 
with galaxies in their sample), to which contributes both the central spheroid and the polar 
disk. We found that NGC4650A is located inside the spread of the data-points and contrary 
to its high luminosity it has metallicity lower than spiral galaxy disks of the same total 
luminosity. If we take into account the total luminosity of the polar disk alone {Mb = —17), 
NGC4650A falls in the region where HII and irregular galaxies are also found, characterized 
by lower luminosity and metallicities with respect to the spiral galaxies. 
For what concerns the chemical abundances in PRGs, onl y few measurements are ava ilable 
in the literature. By using the direct method, very recently iPerez-Montero et al.l (120091 ) have 
derived the chemical abundances of IIZw71, a blue compact dwarf galaxy also catalogued as 
a probable polar ring: consistently with its low luminosity, the metallicity of the brightest 
knots in the ring is lower with respect to that of NGC4650A. Through a similar study, 
Brosch et al.l (120091 ) have derived the chemical abundances for the apparent ring galaxy SDSS 
J075234. 33-1-292049. 8, which has a more similar morphology to that observed in NGC4650A: 
the average value for the oxygen abundance along the polar structure is 12 -|- logO/H = 
8.49 ±0.08. As pointed out by authors, taking into account the ring brightness, such value is 
somewhat lower than that expected by the metallicity-luminosity re lation. A spectroscopi c 
st udy of the peculia r galax y UGC5600, classified as probable PRG ( IWhitmore et al.lll990l ). 
by lShalyapina et al.l (120021 ) reports a quasi-solar metallicity for such object, 12 + log{0 / H) ~ 
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8.8 = 0.9Zq, which is consistent with its high luminosity [Mb 
respect to the values of other PRGs given before. 



-19.4), but it is larger with 



- Metallicity gradient along the polar disk 

The oxygen abundance in the polar disk of NGC4650A as a function of the radius derived 
by empirical and direct methods are shown in Fig. [6] and Fig. [8] respectively: both of them 
show that metallicity remains constant along the polar disk. This suggests that the star 
formation and metal enrichment in the polar structure of NGC4650A is not influenced by 
the stellar evolution of the older central spheroid, w here the last burst of star formation 
occurred between 3 to 5 Gyrs ago (llodice et al.ll2002l ): this turns also to be consistent with 
a later formation of the polar d isk. The absence of any meta llicity gradient is a typical 
behavior found in LSB galaxies ( de Blok fc van der Hulstlll998l ). and also in the polar disk 
galaxy studied by iBrosch et al.l (120091 ). which has a very similar structure to NGC4650A. 
On the contrary, ordinary and oxygen-rich spiral ga laxies show a decreasing abundance with 
increasing radii (see Fig. [6] and iPilyugin et al.ll2006l ). These observed features in spiral disks 
are well explained by the infall mode ls of galaxy formation which predict that they build up 
via accretion by growing inside-out ( iMatteucci fc Francois! Il989t iBoissier fc Prantzod 119991 ) 
and such process generates the observed gradients. 

- Star Formation Rate ( SFR ) 

We have derived th e SFR for the polar disk, from the Ha luminosity using the expression 
given bv lKennicuttI fll998h SFR{MQ/yr) = 7.9 x 10"^^ x L{Ha){erg / s). We found that it is 
almost constant along the disk, within a large scatter of the datapoint. Therefore, given the 
average value of L{Ha) — 3.8 x 10^^ erg/s we have obtained an average SFR ~ 0.04MQ/yr 
and SFR ~ 0.02MQ/yr for the North and South arms respectively. T hese values are quite 



simil ar to the total SFR derived for the HII regions of PRG IIZw71 (jPerez-Montero et al. 



20091 ). which is SFR ~ 0.03 5Mc7)/yr. In order to make a co mparison with SFR estimates for 
a sample of PRGs studied by lReshetnikov fc CombesI (119941 ). we have derived a mean value of 
the Ha flux by integrating within a rectangular aperture of 2.5" x 4.6" along each arm of the 
polar disk and we have obtained an average value of SFR ~ 0. 006 Mp^ / yr /pc^. This value is 



consis tent with the SFR of the HII regions in the PRGs studied by iReshetnikov fc Combes 
Jl994h . where 0.001 < SFR < 0.31. In particular, two PRGs of this sample, UGC7576 
and UGC9796, where the polar structure is very similar to that of NGC4650A (i.e. ex- 
ponential surface brightness profile, very blue colors, knotty appearance, prominent HII 
regions and a large amount of HI gas, all associated to this com po nent and distribu ted as 
a disk in differential rotation; see IReshetnikov &: CombesI Il994j and ICox et al.l 120061 ). have 
SFR ~ 0.005MQ/yr/pc^. As already stressed by lReshetnikov fc CombesI (119941 ). the HII re- 
gions of the PRGs, inc luding also NGC4650A, have SFR similar to late-type spiral galaxies 
( iKennicutt et al.lll987l ). In particular, by comparing the Ha luminosity versus the total B 
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magnitude, the polar disk galaxies, NGC4650A, UGC7576 and UGC9796, are located in the 
same area where Sc-Irr galaxies are also found. 

Taking into account that the po lar disk is very yo ung, since the last burst of star for- 
mation occurred less than 1 Gyr ago (llodice et al.ll2002l ). we check if the present SFR of 0.06 
M^lyr^ and even 2 and 3 times higher (i.e. SFR = 0.12MQ/yr and SFR = O.lSMo/yr), 
can give the estimated metallicity of Z = 0.2Zq and how strongly could incre ase the 
metallicity with time. We used a linearly declining SFR ( iBruzual fc CharlotI l2003l ) ip{t) = 
2M^r~^[l — (t/r)] (typically used for late-type galaxies), to estimate the expected stellar mass 
for the three different values of the SFR (0.06, 0.12 and 0.18 MQ/yr) and three epochs (0.8 
Gyr, 1 Gyr and 2 Gyrs), obtaining stellar masses in the range 4 x lO^ Mfr^ < M^, < 1 x lO-^ ^MtD. 



The stellar mass (M^ ~ 4x lO^M©) in the disk from NIR observations (llodice et al. Il2002h falls 



within this range. Than, by using the mass-metallicity relation derived by iTremonti et al. 
fc004l ). where 12 + log{0 / H) = -1.492 + 1.847/o^f(M^) - 0.08026(/oc/M^)2, we found that 



1.02Zq < Z < IAZq. This shows that the present SFR for the polar disk {SFR = 
OMMQyr-^) is able to increase the metallicity of about 0.2Zq after IGyr. The derived 
values for Z are larger than Z = (0.2 ± 0.002)^0, found by using the element abundances 
(see Sec J3.ip : this differences could be attributed to the accretion of metal-poor gas, as 
discussed in detail in the following sections. 



4.2. Theoretical predictions 



How galaxies acquire their gas is an open issue in the models of galaxy formation and 
recent theoretical works have argued that cold accretion plays a major role. This idea is 



anism in the galaxy formation (IKatz &: White 
Dekel fc Birnboimlbood: lOekel fc Birnboim 2008 



1993I: 


Katz et al. 


1994 




Keres et al. 


2005; 


Bournaud & Elmeereen 


2009). 


Keres et al. 



(I2OO5I ) studied in detail the physics of the cold mode of gas accretion and one feature they 
found is that it is generally directed along filaments, allowing galaxies to draw gas from large 
distances. 



Recently lOcvirk et al.l ( l2008l ). using cosmological simulations, studied the metal enrichment 
of the intergalactic medium as a key ingredient in determining the transition mass from cold 
to hot domin ated accretion. Th e ir me asurements turn to be consistent with the analytical 
prediction of iDekel fc BirnboimI (120061 ) in the low metallicity regime for the cold streams. 
The efficiency of radiative cooling decreases towards low metallicity, and at Z/Zq = 10~^ 
the cooling properties of the gas are those of a primordial mixture. Gas accretion is bi- 
modal both in temperature and in metallicity: the cold accretion mode is associated with a 
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combination of metal-poor filamentary accretion and dense metal-rich satellite galaxy disc 
surrounding, while the hot mode features has strong chemical heterogeneity and a radius- 
dependent metallicity. 

More re c ent si mulations of disk formation in a cosmological context performed by 
Agertz et al.l (120091 ) reve aled that the so calle d chain-galaxies and clump-clusters, found 
only at higher redshifts (lElmegreen et al.l 120071 ). are a natural outcome of early epoch en- 
hanced gas accretion from cold dense streams as well as tidally and ram-pressured stripped 
material from minor mergers and satellites. This freshly accreted cold gas settles into a large 
disk-like systems; simulations show that the interaction region between the new-formed disk 
and the cold st reams can also resul t not aligned with the initial galactic disk: based on a very 
poor statistics, lAgertz et al.l (|2009[ ) suggest that this misalignment might not be typical, and 
it is due to a third cold stream that is perpendicular to the main filament. A more recent 
analysis by iDekel et al.l (120091 ) shows that the accretion of gas along misaligned filaments 
with respect to the disk plane are more common and it leaves traces down to low redshift. 
A n almost polar filame nt can result just as an extreme case of such process and, as suggested 
by lAgertz et al.l ( l2009l ). it could be responsible for the formation of polar disks. 
This scenario not only reproduces the observed morphology and gjlobal rotation of dis ks, but 
also finds a realistic metallicity gradient and SFR of 20MQ/yr. lAgertz et al.l (|2009[ ) found 
solar metallicity for the inner disk, while that in the clump forming region is only ~ I/IOZ© 
due to the accretion of pristine gas in the cold streams mixing with stripped satellite gas. 
Therefore, the studies cited above suggest that the chemical abundance is one of the key 
parameters that can be estimated in a galaxy disk and directly compared with the theoretical 
predictions. This is the main goal of the present work: we studied the chemical abundances 
in the polar disk of NGC4650A in order to check the cold accretion scenario for this object. 



Hydrodynamical simulations performed by iMaccio et al.l ( l2006l ) and iBrook et al.l (120081 ) 
have shown that the formation of a polar disk galaxy can occur naturally in a hierarchi- 
cal universe, where most low-mass galaxies are assembl ed through the accr etion of cold gas 
infalling along a filamentary structures. According to iMaccio et al.l (|2006[ ). the polar disk 
forms from cold gas that fiows along the extended ~ IMpc filament into the virialized dark 
matter halo. The gas streams into the center of the halo on an orbit that is offset from radial 
infall. As it reaches the center, it impacts with gas in the halo of the host galaxy and with 
the warm gas flowing along the opposite fllament. Only the gas accreted perpendicular to 
the major axis of t he potential can survive for more than a few dynamical lifetimes. 
Brook et al.l ( l2008l ) argued that polar disk galaxies are extreme examples of the misalignment 
of angular momentum that occurs during the hierarchical structure formation: an inner disk 
starts forming shortly after the last major merger at z ~ 2. Due to its gas rich nature, 
the galaxy rapidly forms a new disk whose angular momentum is determined by the merger 
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orbital parameters. Later, gas continues to be accreted but in a plane that is almost perpen- 
dicular to the inner disk. At z ~ 0.8 the central galaxy is still forming stars in a disk, while 
the bulk of new star formation is in the highly inclined polar disk. By z ~ 0.5 the inner disk 
has exhausted its gas, while gas continues to fall onto the polar disk. From this point, star 
formation occurs exclusively in the polar disk, that remains stable for at least 3 Gyrs. 
The formation mechanisms described above can self-consistently explain both morphology 
and kinematics of a polar disk galaxy and, in particular, all the observed features (like col- 
ors and color gradients, longevity, spiral arms, HI content and distribution) of the polar 
structure. 



5. Summary and conclusions 



The present study could be considered a step forward both to trace the formation history 
of NGC4650A and to give hints on the mechanisms at work during the building of a disk by 
cold accretion process. 

As mentioned in the SecJH the new kinematic data obtained for the central spheroid 
( llodice et al.ll2006l ). together with the previous studies, set important constraints on the 
possible formation mechanisms for NGC4650 A. In particular, the rn erging scenario is ruled 



out because, according to simulations (e.g. iBournaud et al.l l2005l ). a high mass ratio of 



the two merging galaxies is required to form a massive and extended polar disk as ob- 
served in NGC 4650A: this would convert the intruder into an elliptical- like, not rotationally 
supported, stellar system. This is in contrast with the high maximur n rotation velocit y 
(~ 80 -7- 100 km/s) observed in the outer regions of the central spheroid (jlodice et al.ll2006l ). 



Both the high baryonic mass (star plus gas) in the polar structure and its large extension 
cannot reconcile with a polar ring formed via the gradual disruption of a dwarf satellite 
galaxy (as explained by llodice et al.l 120021 ). Differently, a wide polar ring and/or disk (as 
observed in NGC4650A) may form both around a disk or an elliptical galaxy through the 
tidal accretion of gas stripped from a gas-rich donor, in the case of larg e relative velocity and 



large impact parameter and for a particular orbital configuration (e.g. iBournaud fc Combes 
lioos). Therefore, the two formation scenarios which can be really envisioned in the specific 



(Maccio et al. 


2006; 


Brook et al. 


2008) 



derived the metallicity and SFR for the polar disk in NGC4650A in order to compare them 
with those predicted by different formation scenarios. 



The main results of the present work are (see also Sec. 14. ip : %) the low value of the 
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metallicity derived for the polar disk Z = 0.2Zq, in spite of its high luminosity, Mb = —19.3 
(see Fig. [9]), ii) the lack of any metallicity gradient along the polar disk, which suggests that 
the metal enrichment is not influenced by the stellar evolution of the older central spheroid. 
Hi) the metallicities expected for the present SFR at three different epochs, 1.02Zq < Z < 
IAZq, are higher than those measured from the element abundances and this is consistent 
with a later infall of metal-poor gas. 

In the following we will address how these results reconcile with the predictions by the- 
oretical models (see Sec. 14. 2p and may discriminate between the two formation mechanisms. 



If the polar ring/disk formed by the mass transfer from a gas-rich donor galaxy, the 
accreted material comes only from the outer and more metal-poor regions of the donor: is 
the observed metallicity for the polar component in NGC4650A consistent with the ty pical 
values for the outer disk of a bright spiral galaxy? According to iBresolin et al. koO^ . the 
metallicity of very outer regions of a bright spiral is 0.2Zq < Z < 1.1 Zq: the observed value 
for NGC4650A, Z = 0.2Zq, is close to the lower limit. 

The cold accretion mechanism for disk forma t ion predicts quite lo w metallicity {Z = O.IZq) 
( Dekel &: BirnboimI l2006l . lOcvirk et al.l l2008l . lAgertz et al.l 120091 ) : such value refers to the 
time just after the accretion of a misaligned material, so it can be considered as initial value 
for Z before the subsequent enrichment. How this may reconcile with the observed metallicity 
for NGC4650A? We estimated that the present SFR for the polar disk {SFR = OMMQyr--'^) 
is able to increase the metallicity of about 0.2 after IGv r fsee SecJ4.ip: ta king into account 



that the polar structure is very young, less than IGyr (llodice et al.l 120021 ). an initial value 
of Z = O.IZ0, at the time of polar disk formation, could be consistent with the today's 
observed metallicity. 

This evidence may put some constraints also on the time-scales of the accretion process. 
The issue that need to be addressed is: how could a cosmic flow form a ring/disk only 
in the last Gyr, without forming it before? Given that the average age of 0.8 Gyr refers 
to the last burst of star formation, reasonably the gas accretion along the polar direction 
could have started much earlier and stars formed only recently, once enough gas mass has 
been accumulated in the polar disk. Earlier-on, two possible mechanisr as may be proposed . 
One process that may happen is something similar to that suggested by iMartig et al.l (120091 ) 
for the quenching of star formation in early-type galaxies: given that star formation takes 
place in gravitationally unstable gas disks, the polar structure could have been stable for a 
while and this could have quenched its star formation activity, until the accumulated mass 
gas exceeds a stab ility threshol d and star formation resumes. Alternatively, according to 
the simulations by iBrook et al.l ( l2008l ). the filament could have been there for several Gyrs 
with a relative low inclination, providing gas fuel to the star formation in the host galaxy 
first, about 3 Gyrs ago. Then large-scale tidal fields can let the disk/filament misalignment 
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increase over the time till 80-90 degrees and start forming the polar structure during the last 
1-2 Gyrs. This picture might be consistent not only with the relative young age of the polar 
disk, but also with the estimate of the last burst of star formation in the central spheroid 
dlodice et al]l2002h . 



One more hint for the cold accretion scenario comes from the fact that the metallicity 
expected by the present SFR t urns to be h i gher t han those directly measured by the chemical 
abundances. As suggested by iDalcantonI (120071 ). both infall and outflow of gas can change 
a galaxy's metallicity: in the case of NGC4650A a possible explanations for this differenc e 
could be the infall of pristine gas, as suggested by ( iFinlator et al.l 120071 : Ellison et al.ll2008l ). 



The lack of abundance gradient in the polar disk, as typically observed also in LSB 
galaxies, suggests that the picture of a che mical evolution from inside-out, that well repro- 
duce the observed features in spiral disks (IMatteucci &: Francois! Il989l . iBoissier fc Prantzos 
19991 ). cannot be applied to these systems. In particular, observations for LSB galaxies are 
consistent with a quiescent a nd sporadic chemical evolution , but several explanation exist 
that supports such evidence (Ide Blok fc van der Hulstlll998l ). Among them, one suggestion 
is that the disk is still settling in its final configuration and the star formation is triggered 
by external infall of gas from larger radii: during this process, gas is slowly diffusing inward, 
causing star formation where conditions are favorable. As a consequence, the star formation 
is not self-propagating and the building-up of the disk would not give rise to an abundance 
gradient. The similarities be tween LSB galaxies and the polar disk in NGC4650A, includ- 
ing colors (jlodice et al. 2002) and chemical abundances (this work), together with the very 
young age of the polar disk, the presence of star forrning regions toward s larger radii, the 
warping structure of outer arms ( iGallagher et al.ll2002l . Ilodice et al.ll2002l ). and the constant 
SFR along the disk (this work), suggest that the infall of metal-poor gas, through a similar 
process described above, may reasonably fit all these observational evidences. 

Given that, are there any other observational aspects that can help to disentangle in a 
non ambiguous way the two scenarios? 

One important feature which characterize NGC4650A is the high baryonic (gas plus 
stars) mass in the polar structure: it is about 12 x IO^Mq, which is comparable with or 
even higher than the total luminous mass in the central spheroid (of about 5 x lO^M©). In 
the accretion scenario, the total amount of accreted gas by the early-type object is about 
10% of the gas in the disk donor galaxy, i.e. up to 10^ Mq, so one would expect the 
accreted baryonic mass (stellar + gas) to be a fraction of that in the pre-existing galaxy, 
and not viceversa, as it is observed for NGC4650A. Furthermore, looking at the field around 
NGC4650A, the close lumiri ous spiral NGC4650 may be considered as possible donor galaxy 
(iBournaud &: Combed l2003l ). however the available observations on the HI content for this 
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object show that NGC4650 is expected to be gas-poor (lArnaboldi et aLlll997l : Ivan Driel et aL 
20021 ). So, where the high quantity of HI gas in NGC4650A may come from? If the polar 
disk forms by the cold accretion of gas from filaments there is no limit to the accreted mass. 

Given all the evidences shown above, we can infer that the cold accretion of gas by 
cosmic web filaments could well account for both the low metallicity, the lack of gradient 
and the high HI content in NGC4650A. An independent evidence which seems to support 
such scenario for the formation of polar disks comes from the discovery and study of an 
isolated polar disk galaxy, located in a wall between two voids (jStanonik et al.l 120091 ): the 
large HI mass (at least comparable to the stellar mass of the central galaxy) and the general 
underdensity of the environment can be consistent with the cold flow accretion of gas as 
possible formation mechanism for this object. 

The present work remarks how the use of the chemical analysis can give strong con- 
straints on the galaxy formation, in particular, it has revealed an independent check of the 
cold accretion scenario for the formation of polar disk galaxies. This study also confirmed 
that this class of object needs to be treated differently from the polar ring galaxies, where 
the polar structure is more metal rich (like UGC5600, see Sec. 14. ip a nd a tidal accretion or 
a ma jor merging process can reliable explain the observed properties (iBournaud fc Combes 
20031 ) . Finally, given the similarities between polar disks and late- type disk galaxies, except 
for the different plane with respect to the central spheroid, the two classes of systems could 
share similar formation processes. Therefore, the study of polar disk galaxies assumes an 
important role in the wide framework of disk formation and evolution, in particular, for 
what concern the "rebuilding" of disks through accre tion of gas from cosrnic fila ments, as 
predicted by hierarchical models of galaxy formation (jSteinmetz fc Navarro! 120021 ). 
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Table 1: General properties of NGC4650A 



Parameter 


Value 


Ref. 


Morphological type 


PRG 


NED^l 


R.A. (J2000) 


12h44m49.0s 


NED 


Dec. (J2000) 


-40d42m52s 


NED 


Helio. Radial Velocity 


2880 km/s 


NED 


Redshift 


0.009607 


NED 


Distance 


38 Mpc 




Total (niag) 


14.09 ± 0.21 


NED 


M{HI){Mq) 


8 X 10^ 


Arnaboldi et al. ri997) 


Lb{Lq) 


3 X 10^ 





"NASA/IPAC Extragalactic Database 
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Table 2. Observed and de-reddened emission lines fluxes relative to Hj^ in the 3300 — 6210 

A wavelength range. 



r (arcsec)^ 


OII[3727]/i/;3 


H^[4340]/Hp 


OIII[4363]/i/;3 


OIII[4959]/iJ^ 


OIII[5007]/iJ;3 


slitl 


Fobs 


F'int 


Fobs 


F'int 


Fobs 


Fint 


Fobs 


Fint 


Fobs 


Fint 




±0.05 


±0.06 


±0.013 


±0.015 


±0.004 


±0.005 


±0.02 


±0.02 


±0.03 


±0.03 


-44.1 


5.95 


7.76 






0.058 


0.066 


0.52 


0.51 


1.83 


1.77 


-42.34 


4.32 


5.64 






0.201 


0.227 


0.52 


0.51 


1.83 


1.77 


-40.57 


4.01 


5.22 


0.392 


0.445 


0.182 


0.205 


0.58 


0.57 


1.97 


1.91 


-38.81 


3.69 


4.81 






0.094 


0.107 


0.57 


0.56 


1.67 


1.61 


-37.04 


3.61 


4.71 


0.438 


0.497 


0.096 


0.109 


0.80 


0.78 


1.91 


1.85 


-35.28 


3.39 


4.42 


0.398 


0.451 


0.095 


0.108 


0.80 


0.78 


1.82 


1.76 


-33.52 


3.42 


4.45 


0.329 


0.373 


0.064 


0.073 


0.79 


0.78 


1.87 


1.81 


-31.75 


3.38 


4.41 


0.377 


0.428 


0.043 


0.048 


0.51 


0.50 


1.60 


1.55 


-29.99 


3.47 


4.52 


0.423 


0.480 


0.038 


0.043 


0.54 


0.52 


1.58 


1.53 


-28.22 


3.99 


5.21 


0.464 


0.526 


0.188 


0.212 


0.60 


0.59 


1.76 


1.71 


-i ( .04 


4.06 


5.29 


0.400 


0.454 


0.105 


0.119 


1.05 


1.03 


3.11 


3.01 


-14.11 


3.12 


4.06 


0.376 


0.427 






0.50 


0.50 


1.36 


1.31 


-O.OZo 


4.35 


5.67 










0.81 


0.80 


2.51 


2.43 


26.46 






0.348 


0.395 














28.22 






0.540 


0.613 














37.04 


3.39 


4.42 










1.79 


1.75 


5.33 


5.16 


38.81 


3.94 


5.13 










1.92 


1.88 


6.85 


6.63 


45.86 


2.41 


3.14 










0.24 


0.24 


0.51 


0.49 


47.63 


4.33 


5.64 






0.030 


0.034 


0.21 


0.20 


0.37 


0.36 


49.39 


3.70 


4.82 










0.20 


0.20 


0.33 


0.32 


slit 2 






















-49.49 






0.598 


0.678 














-33.52 


1.97 


2.57 


0.555 


0.630 






0.19 


0.19 


0.43 


0.42 


-31.75 


2.24 


2.91 










0.17 


0.17 


0.54 


0.53 


-3.528 


2.54 


3.31 










0.88 


0.86 


2.89 


2.80 


5.292 


3.59 


4.69 


0.233 


0.264 






0.80 


0.78 


2.69 


2.60 


8.82 


3.28 


4.28 










0.43 


0.42 


1.27 


1.23 
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Tablc 2 — Continued 



r (arcsec)^ Oil [3727] /i/;3 H^[4340]/H^ OIII[4363]/i/;3 OIII[4959]/i/^ OIII[5007]/i//j 



17.64 


3.15 


4.11 


0.594 


0.674 






0.92 


0.90 


2.15 


2.08 


21.17 


2.22 


2.90 










0.88 


0.86 


2.56 


2.48 


22.93 


3.27 


4.26 


0.386 


0.438 






0.51 


0.50 


1.58 


1.53 


24.7 


2.44 


3.19 


0.387 


0.439 






0.93 


0.91 


2.71 


2.62 


26.46 


3.08 


4.01 


0.442 


0.501 






0.83 


0.82 


2.51 


2.43 


28.22 


4.13 


5.38 


0.444 


0.504 






0.75 


0.73 


2.27 


2.19 


29.99 


3.32 


4.33 


0.472 


0.535 






0.64 


0.62 


1.91 


1.85 


31.75 


5.96 


7.77 


0.423 


0.480 






0.83 


0.81 


2.22 


2.15 


33.51 






0.540 


0.613 














44.1 


2.57 


3.35 










0.88 


0.86 


2.84 


2.75 


45.86 


4.14 


5.40 


0.427 


0.484 


0.129 


0.146 


0.09 


0.09 


0.43 


0.42 


47.63 


2.02 


2.64 






0.537 


0.606 


0.73 


0.71 


2.27 


2.19 


49.39 






0.384 


0.436 














56.45 


3.46 


4.51 










1.04 


1.02 


3.72 


3.60 


58.21 


3.93 


5.13 






0.162 


0.183 


1.15 


1.13 


4.18 


4.05 


59.98 


4.92 


6.41 










1.68 


1.65 


6.63 


6.42 



^The negative values are for the northern regions of the spectra 
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Table 3. Observed and de-reddened emission lines fluxes relative to i^^ in the 

5600 - llOOOA wavelength range. 



r (arcsec)* 


i/„[6563]/V 


Sll[6717]/Hp 


Sll[6731]/Hfi 


SIII[9069]/i/^ 


SIII[9532]/i7;3 


sht 1 


Fobs 


F'int 


Fobs 


Fint 


Fobs 


Fint 


Fobs 


Fint 


Fobs 


Fint 




±0.02 


±0.03 


±0.015 


±0.011 


±0.016 


±0.012 


±0.08 


±0.06 


±0.15 


±0.11 


-38.81 


















0.05 


0.03 


-37.04 


5.46 


4.25 


















-35.28 


3.85 


2.99 














0.017 


0.010 


-33.52 


2.44 


1.896 














0.04 


0.03 


-31.75 


0.930 


0.724 


















-29.99 


0.412 


0.321 


















-26.46 


















0.05 


0.03 


-17.64 


















0.03 


0.019 


-15.88 


















0.06 


0.03 


-10.58 






0.199 


0.153 


0.632 


0.483 










-8.82 






0.426 


0.329 


0.577 


0.441 


0.052 


0.0320 






-7.06 






0.84 


0.646 


0.935 


0.714 


0.0293 


0.0178 






-5.29 














0.041 


0.0252 






-3.53 














0.68 


0.41 






-1.76 














0.266 


0.162 






0.0 






1.39 


1.07 


1.31 


1.00 










1.76 






1.48 


1.14 


0.609 


0.465 


0.33 


0.202 






3.53 






0.686 


0.529 


0.460 


0.351 










5.29 






0.136 


0.105 


0.422 


0.322 










7.06 






0.0142 


0.0109 


0.377 


0.288 


0.041 


0.0249 


0.023 


0.014 


29.99 


















0.014 


0.008 


31.75 


















0.022 


0.013 


40.57 






0.0340 


0.0262 


0.174 


0.133 










51.16 






0.465 


0.358 


0.280 


0.214 


0.33 


0.202 


1.0 


0.6 



sht 2 



-45.86 0.11 0.06 

-38.81 1.59 1.24 ••■ •■• 
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Tablc 3 — Continued 



r (arcsec)^ i/„[6563]/i/^^ SII[6717]/if^ SII[6731]/i/^ SIII[9069]/i/^ SIII[9532]/i/^ 



-37.04 


5.48 


4.26 


















-33.52 


















0.13 


0.08 


-22.93 


0.874 


0.680 


















-17.64 














0.240 








-5.29 






0.0536 


0.0413 


0.271 


0.207 


■ • • 




0.22 


0.1319 


-3.53 






0.090 


0.070 


0.0551 


0.0421 






. . . 


. . . 


-1.76 






0.962 


0.741 


0.0023 


0.00177 






0.17 


0.10 


0.0 






0.0691 


0.0533 


0.113 


0.087 










1.76 


















0.13 


0.08 


3.53 














0.55 


0.33 


0.3 


0.17 


5.29 














0.36 


0.218 






8.82 






0.583 


0.449 


0.87 


0.658 


1.59 


0.97 


2 


1.3 


10.58 














0.67 


0.41 


0.3 


0.20 


19.4 


















0.18 


0.10 


21.17 














0.54 


0.33 






24.7 


1.63 


1.26 


0.159 


0.123 


1.15 


0.88 






0.18 


0.11 


26.46 


0.290 


0.226 


















28.22 


0.0620 


0.0482 














0.15 


0.09 


29.99 


0.0938 


0.0730 










0.51 


0.311 






31.75 


3.82 


2.97 


0.118 


0.091 


0.291 


0.222 






0.09 


0.05 


33.52 


5.47 


4.25 


0.204 


0.157 


0.216 


0.165 






0.23 


0.13 


35.28 






0.096 


0.074 


0.237 


0.181 










38.81 














0.43 


0.259 






42.34 






















56.45 






0.091 


0.070 


0.230 


0.175 










58.21 


















0.22 


0.13 


59.98 






0.077 


0.0587 


0.134 


0.102 


0.37 


0.226 
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'^The negative values are for the northern regions of the spectra 
^Ha flux is corrected for corresponding [A^//]A6548 Hux/H^ 
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Table 4: Oxygen and Sulphur abundance parameters and oxygen abundances for the polar 
disk in NGC4650A ^^^^^^^^^^^^^^^^^^^^^ 

r (arcsec) /?23 (±0.05) 12 + 1(>g(0 / H)p (±0.2) 



Slit 1 






-44.1 


10.05 


7.8 


-42.34 


7.91 


8.0 


-40.57 


7.70 


8.1 


-38.81 


6.98 


8.1 


-37.04 


7.35 


8.1 


-35.28 


6.97 


8.2 


-33.52 


7.04 


8.2 


-31.75 


6.45 


8.2 


-29.99 


6.57 


8.2 


-28.22 


7.50 


8.1 


-17.64 


9.33 


8.0 


-14.11 


5.87 


8.2 


-3.528 


8.90 


8.0 


37.04 


11.33 


8.0 


38.81 


13.64 


7.8 


45.86 


3.86 


8.4 


47.63 


6.20 


8.0 


49.39 


5.34 


8.1 


Slit 2 






-33.52 


3.18 


8.5 


-31.75 


3.61 


8.4 


-3.528 


6.97 


8.3 


5.292 


8.08 


8.1 


8.82 


5.93 


8.2 


1 7 64 


7 1 


8 2 


21.17 


6.24 


8.4 


22.93 


6.28 


8.2 


24.7 


6.73 


8.3 


26.46 


7.26 


8.2 


28.22 


8.30 


8.0 


29.99 


6.81 


8.2 


31.75 


10.74 


7.7 


44.1 


6.96 


8.3 


45.86 


5.91 


8.1 


47.63 


5.55 


8.4 


56.45 


9.13 


8.1 



58.21 10.30 8.0 
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Table 5: Oxygen abundances for the polar disk in NGC4650A derived by the Tg method 
r(arcsec) 12 + log{0 / H)t^^^ 12 + log{0 / H)t^S^Q.I) 



Slit 1 

-44.10 7.83 ± 0.58 8.6 

-42.34 7.71 ±0.52 8.5 

-40.57 7.69 ± 0.49 8.5 

-38.81 7.65 ±0.50 8.4 

-37.04 7.67 ±0.46 8.4 

-35.28 7.64 ±0.46 8.4 

-33.52 7.65 ± 0.45 8.4 

-31.75 7.62 ±0.49 8.4 

-29.99 7.63 ± 0.49 8.4 

-28.22 7.69 ±0.50 8.5 
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Fig. 1. — Optical image of NGC4650A with superimposed the sUts used to acquire data 
analyzed in this work. 



- 35 - 




Fig. 2. — 1-D spectrum of the standard star LTT4816, used to flux-calibrate the spectra. 
The black line represent the spectrum of the star acquired in the blue wavelength range, 
while the red line is the same star acquired in the red range. 
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Fig. 3. — Top panel - Spectrum of NGC4650A in the blue wavelength range, corresponding 
to the region marked as A in the Fig. [H Middle panel - Spectrum of the region B in the 
Fig. [H Bottom panel - Spectrum of the region C in the Fig. [1], where are clearly visible the 
absorbing features of the line and Hp. 
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Fig. 4. — Top panel - Spectra of NGC4650A in the red wavelength range, corresponding to 
the region marked as A in the Fig. [H We divided the spectra in Al and A2 in order to 
obtain a better visuahzation of the hnes. Bottom panel - Spectrum of the region B (Bl and 
B2 as explained below) in the Fig. [H 
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Fig. 5.— Oxy gen abundance, obtained by using the empirical calibration introduced by 
PilyuginI (l200l[ l. versus oxygen abundance parameter R23, (see Sec. 13. ip . 
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Fig. 6. — Oxygen abundance derived with empirical methods proposed by iPilyuginl (j200l[ ) 
versus radius. The values in Ta b. HI were binned each 5 arcseconds. The superimposed lines 
are the linear best fit derived by lPilyugin et al.l (120061 ): the red line represents the best fit to 
the abundance of oxygen-rich spirals, while the blue line is those related to ordinary spirals. 
The black line is the best fits obtained for NGC4650A. 
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Fig. 7. — Electron temperature along the polar disk of NGC4650A. The filled circles corre- 
spond to the observed T[OIII], the starred points correspond to the T[OII] estimated by the 
photoionization model (see Sec. 13.21 for details), the open square corresponds to the average 
value of T[OII] directly estimated by the ratio Ron = O/J [3727]/ [7325]. In the bottom left 
side of the plot the mean error is shown. 
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Fig. 8. — Oxygen abundance derived with direct method. Squares are the values derived 
by using the observed T[OII] estimates; stars are those derived by using the T[OII] from 
photoionization models. See Sec. 13.21 for details. The errors on oxygen abundance is 0.5. 
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Fig. 9. — Oxygen abundance vs absolute blue magnitude for CNELGs (yellow filled circles), 
ELGs (cyan open circles), four field galaxies with emission lines (filled black circles), nearby 
dwarf irregulars (open triangles), local spiral galaxies (open circl es), local HII galaxies (ope n 
squares), NGC4650A (star) and the polar disk galaxy IIZw71 (jPerez-Montero et al.ll2009l ). 
The dashed line indicates the solar oxygen abundance. The arrow indicates the shift of the 
value of the oxygen abundance if we use the direct methods to evaluate it. The total B-band 
magnitude for NGC4650A (Mr = —19.3 ) has been evaluated by using the same value of 
Hq used by iKobulnicky fc ZaritskyI (119991 ) in order to compare NGC4650A with galaxies in 
their sample 



